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Abstract 
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We study spontaneous emission of an atom near a nanofiber with two fiber-Bragg-grating (FBG) 
mirrors. We show that the coupling between the atom and the guided modes of the nanofiber can 

(N 

be significantly enhanced by the FBG cavity even when the cavity finesse is moderate. We find 

Si 

that, when the fiber radius is 200 nm and the cavity finesse is about 30, up to 94% of spontaneous 
emission from the atom can be channeled into the guided modes in the over damp ed-cavity regime. 
We show numerically and analytically that vacuum Rabi oscillations and strong coupling can occur 
in the FBG cavity even when the cavity finesse is moderate (about 30) and the cavity length is 
large (on the order of 10 cm to 1 m), unlike the case of planar and curved Fabry-Perot cavities. 
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I. INTRODUCTION 



Coupling of light to subwavelength structures and its control pose one of the greatest 
challenges of recent research [l, 2, 3, 4, 5, 6, 
at microwave frequencies has been observed 

enables strong coherent coupling_between individual emitters and guided plasmon excitations 
in conducting nanostructures 



7\ . Strong coupling in a superconducting circuit 
[2]. Chang et al. have proposed a technique that 



In the case of dielectric waveguides, it has been shown 
that a significant fraction (up to 28%) of emission from a single atom can be channeled into 

a — a a a a. — ^ « mjm ■» „ * . ^ b ^ 



studied in the context of a two-level atom 



atom with a hyperfine structure of energy levels [4j. 



as well realistic multilevel cesium 
The parameters for the decay of 

ated 



Eopulations and cross-level coherences of an atom near a nanofiber have been calcu 
], 0]. The cooperation of distant atoms via a nanofiber has been discussed 11, 12|. It 
has been shown that, due to guided modes, a substantial cooperation can survive large 
interatomic distances ll| . and a linear array of distant atoms can significantly enhance the 
rate of spontaneous emission and the efficiency of channeling of emission into the nanofiber 
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Optical cavities are often employed to increase the interaction between atoms and photons 
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Various cavity quantum 



electro dynamic effects have been studied [13]. There have been spectacular recent successes 



brought by the merging of optical cavity systems with ultracold neutral atoms 



23. 



1 

241 



19 

25 



20 



26 



21 



23, 



0, 



22 1 as well as with electromagnetically induced transparency physics 
28l . l29j . It is natural to expect that the use of a cavity can substantially 



enhance the channeling of emission from an atom into a nanostructure. It is desirable to 
combine the cavity technique with the nanofiber technique to obtain a hybrid system, where 
the interaction is enhanced by the transverse confinement of the field in the fiber cross- 
section plane as well as the longitudinal confinement of the field between the mirrors. Such 
a system has been studied recently in the context of intracavity electromagnetically induced 
transparency [3^]. It has been shown that the presence of a fiber-Bragg-grating (FBG) 
cavity with a large length (on the order of 10 cm) and a moderate finesse (about 30) can 
significantly enhance the group delay of the guided probe field j^ . 

In this paper, we study spontaneous emission of an atom near a nanofiber with two 
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FBG mirrors. We find that the coupling between the atom and the guided modes can be 
significantly enhanced by the FBG cavity even when the cavity finesse is moderate. We 
show numerically and analytically that vacuum Rabi oscillations and strong coupling can 
occur in the FBG cavity even when the cavity finesse is moderate (about 30) and the cavity 
length is large (on the order of 10 cm to 1 m). 



Before we proceed, we note that there has been a 



Bragg gratings over the past two decades 



31. 
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arge body of work involving fiber 
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361 ] . With careful control of the 



grating writing process and appropriate choice of glass material, a FBG resonator can have 
a finesse of well over 1000 and a linewidth of a few MHz [36]. 

The paper is organized as follows. In Sec. [TT]we describe the model of a nanofiber with two 
FBG mirrors. In Sec. IHII we derive a basic equation for spontaneous emission of an atom in 
the model. In Sec. IIVI we study spontaneous emission of the atom in the overdamped-cavity 
regime. In Sec. |V]we derive a delay-differential equation for spontaneous emission and study 
it numerically. In Sec. IVII we approximate the delay-differential equation under the single- 
mode cavity condition and analyze the atomic decay in various cases. Our conclusions are 
given in Sec. IVHI 

II. MODEL 

We consider spontaneous emission of a two-level atom in the vicinity of a nanofiber with 
two FBG mirrors (see Fig. [I]). The field in the guided modes of the nanofiber is reflected 
back and forth between the FBG mirrors. Such a system is a nanofiber-based cavity. The 
nanofiber has a cylindrical silica core of radius a and of refractive index n\ = 1.45 and an 
infinite vacuum clad of refractive index n 2 = 1. In view of the very low losses of silica in the 
wavelength range of interest, we neglect material absorption. We also neglect the effects of 
the surface-induced potential, the surface roughness, and the phonon heating on the atom. 
We use the cylindrical coordinates (r, (p, z), with z being the axis of the fiber. 

In the presence of the fiber, the electromagnetic field can be decomposed into guided and 
radiation modes 3t|]. In order to describe the field in a quantum mechanical treatment, we 
follow the continuum field quantization procedures presented in j^]. First, we temporally 
neglect the presence of the FBG mirrors. Regarding the guided modes, we assume that the 
single-mode condition is satisfied for a finite bandwidth around the atomic transition 
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FIG. 1. (Color online) Spontaneous emission of an atom in the vicinity of a nanofiber with two 
fiber-Bragg-grating mirrors. 

frequency u> . We label each fundamental guided mode HE n with a frequency u in this 
bandwidth by an index \i = (u,f,l), where f = +,— denotes the forward or backward 
propagation direction and I = +, — denotes the counterclockwise or clockwise rotation of 
polarization. In the interaction picture, the quantum expression for the electric positive- 
frequency component of the field in the cavity-free guided modes is 4] 



Here = e^) (r, ip) is the profile function of the guided mode \i in the classical problem, 
a M is the corresponding photon annihilation operator, ^ = J °° dw is the summation 
over the guided modes, (3 is the longitudinal propagation constant, and (3' is the derivative 
of f3 with respect to uj. The constant j3 is determined by the fiber eigenvalue equation \^\ . 
The operators a M and satisfy the continuous-mode bosonic commutation rules [a M , a' ,] = 



8(u - u')5ff>5u>. 



{ see also Refs. {4. l37l|). 

Regarding the radiation modes, the longitudinal propagation constant (3 for each fre- 
quency lu can vary continuously, from —k to k, with k = u/c being the wave number. We 
label each radiation mode by an index v = (u,/3, m, I), where m = 0, ±1,±2, ... is the 
mode order and I — +, — is the mode polarization. In the interaction picture, the quantum 
expression for the electric positive-frequency component of the field in the radiation 
modes is jj] 

Here e^) = (r, ip) is the profile function of the radiation mode v in the classical problem, 
a„ is the corresponding photon annihilation operator, and = ^2 ml J °° duo f* k d(3 is the 
summation over the radiation modes. The operators a u and a\ satisfy the continuous-mode 



he explicit expression for the mode function is given in Appendix 
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bosonic commutation rules [a u , a],,] = 5(cu — u')6(P — (3')5 mm >8u>. The explicit expression for 
the mode function is given in Appendix iBl (see also Refs. {4. \s7\). 

Next, we take into account the effect of the FBG mirrors on the mode functions. We 
assume that the two FBG mirrors are identical, having the same complex reflection and 
transmission coefficients R and T, respectively, for the guided modes in a broad bandwidth 
around the atomic transition frequency ujq. In general, we have |i?| 2 + |T| 2 < 1, where the 
equality (inequality) occurs for lossless (lossy) gratings. Without loss of essential physics, we 
assume that the gratings are lossless, that is, |i?| 2 + |T| 2 = 1. Let the mirrors be separated 
by a distance L, from the point z = —L/2 to the point z = L/2. The mode functions of the 
guided modes are modified by the presence of the mirrors. The forms of the cavity- modified 



mode functions are obtained, as usual in the Fabry-Perot theory, 



series resulting from the multiple reflections by the mirrors 



39 
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40 



summing the geometric 



41] . Inside the cavity, 



the mode functions of the cavity-modified guided modes are given by 

s(w,+,0 = r fo+,0 1 + c (u-J) IIte 

1 - R 2 e 2i ? L 1 - R?e 2i < 3L ' 

1 - R 2 e 2i f 3L 1 - R 2 e 2i P L ' 



(3) 



and, hence, the electric positive-frequency component of the field in the cavity-modified 
guided modes is 

^cavgyd = * E ^ ^ ^ ' ' W 

We assume that the FBG mirrors do not reflect the radiation modes. This assumption is 
reasonable in the case where the distance L between the FBG mirrors is large as compared 
to the fiber radius a and to the wavelength Ao = 2ir/k , with k Q = ujq/c being the wave 
number of the atomic transition. With this assumption, the mode functions of the radiation 
modes are unchanged by the presence of the FBG mirrors. Inside the cavity, the electric 
positive-frequency component of the total field is given by 

E {+) = EW yd + E« (5) 

We emphasize that the FBG cavity described above confines only the guided modes, 
whose wave vectors are aligned along the fiber axis direction z. The radiation modes are not 
confined by the FBG cavity. In this sense, the physics of the FBG cavity is similar to that 



of one- dimensional cavities 4l|, |42j, and is different from that of planar Fabry-Perot cavities 



13 



39 



40 



431 ]. where off-axis modes reduce the quantum electrodynamic (QED) effect of 



40 



431 ] . We also note that the guided field 



the cavity on spontaneous emission of the atom 
in the FBG cavity is confined not only in the axial direction between the mirrors but also in 
the fiber cross-section plane. In this sense, the physics of the FBG cavity is similar to that 



of curved 
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abry-Perot cavities, which are often used in experiments on cavity QED effects 
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2l| . An advantage of a FBG cavity based on a nanofiber is 



that the field in the guided modes can be confined to a small cross-section area whose size is 
comparable to the light wavelength For example, for a nanofiber with radius of 200 nm, 
the effective mode area A e g = (J |e^| 2 G?r) 2 / J |e^| 4 c/r of the fundamental guided modes 
with the wavelength A = 852 nm is found to be A e $ = 0.65 /im 2 . The corresponding mode 
radius is found to be r e ff = a/ A c r / tx = 454 nm, which is much smaller than the typical values 
of 15 to 30 urn for the waists of the cavity modes used in the experiments on cavity QED 



effects 
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19 



2d, 



211 ] . Another advantage of the nanofiber-based cavity 



is that the cavity guided field can be transmitted over long distances for the communication 
purposes. 

We now describe the interaction between the atom and the field. Let \a) and \b) be the 
upper and lower states of the atom, respectively. The operators o = \b)(a\ and = \a)(b\ 
describe the downward and upward transitions of the atom, respectively. In the dipole 
and rotating-wave approximations and in the interaction picture, the Hamiltonian for the 
atom-field interaction is 



ih G u a ] a u e- i{uj - ua)t + H.c, 



(6) 



where the coefficients and G v characterize the coupling of the atom with the cavity- 
modified guided modes fi = (ui,f,l) and the radiation modes v = (u, f3,m,l), respectively 
Their explicit expressions are 



G, 



G h 



a / — — — rd-e^(r,^)le^ +m ^. 
V 4ire n L J 



(7a) 
(7b) 
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Here d = (a|d|6) is the matrix element of the electric dipole moment of the atom, and r, tp, 
and z are the cylindrical coordinates of the position of the atom. 

III. BASIC EQUATION FOR SPONTANEOUS EMISSION 

We assume that the atom is initially excited and the field is initially in the vacuum state. 
The wave function of the combined atom-field system at an arbitrary time t can be written 

as 

M = C a \a; 0) +J2 C ^\b; U + J2 C »^ ^ 

fx v 

Here C a is the probability amplitude for the atom to remain in the upper state \a), and 
and Cbu are the probability amplitudes for the atom to move to the lower state \b), emitting 
a photon into a guided mode /x and a radiation mode u, respectively. In the interaction 
picture, the Schrodinger equation = H- 1Ilt \ip) yields the following equations for the 

probability amplitudes: 

Ca = ~Y, G^^Cb, - G.e-^-^Cbu (9) 

H v 

and 

Cbu = Gle^-^Ca. (10) 
We integrate Eqs. ffTUj) and substitute the results into Eq. (jUJ). Then, we obtain 
C a(t) = ~ Y\ |G M | 2 [ e-^-^ T C a (t - r)dr 

-Y,\GA 2 I e~ l ^-^ T C a (t-r)dr. (11) 

Since the radiation modes are not confined by the FBG cavity, the interaction between 
the atom and the field in the radiation modes is weak. In addition, the mode functions 
e>( u ' of the radiation modes are smooth with respect to the mode frequencies. Therefore, 
we can apply the Born-Markov approximation to the contribution of the radiation modes, 
that is, to the terms associated with the second integral on the right side of Eq. (TTTT) . 
In this approximation, we replace C a (t — r) by C a (t) and take it out from the integral. 
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With the assumption that the observation time t is much larger than the atomic oscillation 
period 2tt/ujo, we extend the upper integration limit t to +00. Furthermore, we neglect 
the imaginary part of the result of the integration, which describes the contribution of the 
radiation modes to the Lamb shift of the atomic transition frequency. Then, we obtain 

Ca(t)=[C a (t)] gyd -^fC a (t), (12) 

where the term 

[^gyd = ~ E l^l 2 J* e-^-^C a (t - r)dr (13) 
describes spontaneous emission into guided modes and the coefficient 

7rad = 27r^|G i ,| 2 <5(cu-^o) (14) 

V 

is the rate of spontaneous emission into radiation modes. 

In terms of the mode functions e^) of the radiation modes, expression (fT4|) for 7 rac j can 



be rewritten as 



tt-t E r^|d-e^) (r ^)| 2 . (15 ) 



7rad 

ml 



The rate Trad of spontaneous emission into radiation modes has been calculated and studied 



in Refs. 14 



111- 



Since the mode functions of the guided modes are modified by the FBG mirrors, they 
may contain narrow resonances. Therefore, we need to perform a special treatment for the 
contributions from the guided modes. 

We introduce the notation V = V z and V±\ = t(V x ± iV y )/y2 for the spherical compo- 
nents of an arbitrary vector V, and the notation uo = z and u-ti = =p(x ± iy)/\/2 for the 
spherical basis vectors. Without loss of essential physics, we assume that only one spherical 
component d q = d of the dipole vector d, where q = —1, 0, or 1, is nonzero. 

We use Eq. ([3]) to calculate Eq. f!7a|) and then insert the result into Eq. ffTSl . We make 
two approximations. One is to allow the frequency u to be negative for the convenience of 
calculation. The other is that the guided-mode functions = e^'^ and the factor u/3' 
are estimated at the atomic transition frequency uj . These approximations are valid because 
the oscillations described by the exponential factor e -*( w - w o)T j n (J13J) are generally very 
fast except for a small region where the mode frequency uj is close to the atomic transition 
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frequency ujq. As a result, we obtain the following equation for the contribution of the 
cavity- modified guided modes to the atomic decay: 



x C a (t - r)dr. (16) 



Here 

Tgyd = ^f- £ l d ■ « (U0, °r = ^(1-/ + (17) 

2e hv g 1 1 e ftu g 

is the rate of spontaneous emission into guided modes in the absence of the FBG mirrors 
and 

1 + |#| 2 + 2|,R| cos<5 cos(2/fe) 

is the cavity impact (enhancement /inhibition) factor. 

In Eq. ffT7]) . we have introduced the notation |eo| = |ei aJ0 ' +,+ ' ) | and |e±i| = (|er aJ °' + ' + ' ) | =F 
\e^°' + \) / \/2 for the magnitudes of the spherical components of the resonant guided mode 
functions. We have also introduced the notation v g = 1//3'(uj ) for the group velocity of the 



resonant guided field. The cavity-free ratejy. 
has been calculated and studied in Refs. 
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spontaneous emission into guided modes 



In Eq. (1181) . we have introduced the notation 

$ = PL + <j) R + (1 + q)n (19) 

for the shift of the phase of the parallel-to-dipole component of the guided field per cavity 
crossing with a single reflection. Here 4>r is the phase of the complex reflection coefficient R, 
that is, R = \R\e lr ^ R . Depending on the phase shift per cavity crossing $ of the cavity guided 
field and the axial position z of the atom, the cavity impact factor G(u) can be larger or 
smaller than one, indicating enhancement or inhibition, respectively, of spontaneous emission 
into guided modes. Such enhancement and inhibition of spontaneous emission are the Purcell 



effect (45|, which has been studied widely in literature [13]. 

Equation ffl2|) with the term [Ca(£)] gyd given by Eq. ffl6|) is the basic equation for 
spontaneous emission of the atom in the model. We will use this equation to study the 
emission of the atom in different regimes. 
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IV. EXPONENTIAL DECAY IN THE OVERDAMPED-CAVITY REGIME 



We consider the case where the interaction between the atom and the cavity field is weak. 
We assume that the cavity resonance width k is much larger than the characteristic atomic 
decay rate T, that is, k ^> V . In addition, we assume that the observation time t is much 
longer than the atomic oscillation period 2ir/uo, the cavity crossing time tl = L/v g , and the 
cavity damping time k -1 , but is much shorter than the atomic decay time that is, we 
have u t ^> 2ir, t ^> tl, Kt 3> 1, and Tt <C 1. Under these conditions, the Fermi golden rule, 



which is based on the Born-Markov approximation, is valid [46j . We apply the Born-Markov 
approximation to Eq. (1161) for the contribution of the cavity-modified guided modes to the 
atomic decay. In this approximation, we replace C a (t — r) by C a (t). Furthermore, we extend 
the upper integration limit t to +oo. Then, we obtain 

[Ca(t)] gyd = -^C a (t). (20) 

Here 

7cavgyd = Tgyd^o (21) 

is the cavity-modified rate of spontaneous emission into guided modes, with 



Go = GM = - , Vl^r-ri ( 22 ) 



1 + 1 


R 


2 + 2| J R|cos$ 


cos(2/3 ^) 


1 - 


R 


2 


+ ±\R 


l 2 ( 


1 - 


R 


2 ) _1 sin 2 $ 



being the resonant cavity impact factor. In Eq. (1221) . we have introduced the notation 
Po = /3(uj ) for the propagation constant of the resonant guided light and the notation 

$o = $M = /3 L + (f) R + (l+q)7r (23) 

for the resonant-light phase shift per cavity crossing with a single reflection. In deriving 
Eq. (|20|) . we have neglected the contribution of the guided modes to the Lamb shift of the 
atomic transition frequency. Note that expression ff22l) for the resonant cavity impact 



Go is in agreement with the corresponding results for one-dimensional cavities 41 



421. 



'actor 



We insert Eq. ( l20i) into Eq. ( fl2i) . Then, we obtain the exponential-decay equation 

Cait) = ~C a (t), (24) 

with the total atomic decay rate 

T = 7cavgyd + 7rad = lgydG + 7rad- (25) 
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Thus, in the overdamped-cavity regime, the spontaneous emission of the atom is an expo- 
nential decay process. Note that the cavity impact factor Go and, consequently, the rates 
7cavgyd and T depend on the mirror reflection coefficient R. They oscillate with varying z. 
They also oscillate with varying cavity length L through their dependences on the phase 
shift per cavity crossing $ - 

The cavity resonance condition is $o — miT , where m is an integer number. Under this 
resonance condition, Eq. ff22l) for the resonant cavity impact factor Go reduces to 



1 + 


R 


2 + 2\R\ 


cos(2j3 z + rmr) 






1 




\R 


2 



G = ^ ' , D , ™ (26) 



The maximal value 

G max = max(Go) = \ + | p| (27) 
1 — | Jx | 

of the factor Go describes the maximal enhancement of spontaneous emission into guided 
modes. The minimal value 

G min = min(Go) = ^j— f| (28) 
1 + \K\ 

of the factor Go describes the maximal inhibition of spontaneous emission into guided modes. 
It is interesting to note that, under the resonance and overdamped-cavity conditions, the 
maximal enhancement factor G max and the maximal inhibition factor G m ; n do not depend 
on the cavity length L. They depend only on the mirror reflection coefficient R. 

The above results are different from the general results for planar Fabry-Perot cavities 



41 



40, 



43J]. However, they are in agreement with the results for one-dimensional cavities 



, 4^| and also with the results for very narrow planar Fabry-Perot cavities |39j. The reason 
is that the FBG cavity reflects only the fiber guided modes, which propagate along the fiber 
axis, and is therefore similar to one-dimensional cavities. It is known that the enhancement 
factor for a one- dimensional cavity is, in general, larger than that for a corresponding planar 
Fabry-Perot cavity [40]. Therefore, we expect that the FBG cavity can substantially enhance 
the rate of spontaneous emission into guided modes even when the finesse of the FBG cavity 
is moderate. Indeed, for the mirror reflectivity |i?| 2 = 0.8 or 0.9, which correspond to the 
finesse F = 7r|i?|/(l — |i?| 2 ) = 14 or 30, respectively, we obtain the enhancement factor 
G max = 18 or 38, respectively. Such values of the enhancement factor are rather significant. 
For comparison, we note that the maximum enhancement factor for a planar Fabry-Perot 
microcavity with |i?| 2 = 0.9 and L = Xq/2 is just about 3 39j, |4Qj. 
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FIG. 2. Dependences of the total spontaneous emission rate T (solid lines) and its two components 
Tcavgyd (dashed lines) and 7 ra d (dotted lines) on (a) the FBG mirror reflectivity \R\ 2 and (b) the 
fiber radius a in the case where the dipole of the atom is oriented along one of the spherical basis 
vectors u±i. In (a), the fiber radius is a = 200 nm. In (b), the reflectivity of the FBG mirrors 
is \R\ 2 = 0.9. In both cases, the atom is located on the fiber surface (r = a) and at the cavity 
center (z = 0). The length of the cavity is such as the phase shift per cavity crossing <I>o is an even 
multiple of tt. The wavelength of the atomic transition is Ao = 852 nm. The refractive indices of 
the fiber and the vacuum clad are n\ = 1.45 and ri2 = 1, respectively. The rates are normalized to 
the free-space decay rate 70. 

The total spontaneous emission rate Y and its components 7 cav gyd and 7 ra d depend on 
the FBG mirror reflectivity |i?| 2 and the fiber radius a. In Fig. [2] we plot T, 7 cav gyd, and 
7 ra d as functions of |i?| 2 and a in the case where the dipole of the atom is oriented along 
one of the spherical basis vectors u±i. The atom is located on the fiber surface and at 
the cavity center. The length of the cavity is such that the phase shift per cavity crossing 
$0 is an even multiple of ir, that is, an even-order resonance is produced. Under this 
resonance condition, the center of the cavity corresponds to an antinode of the parallel-to- 
dipole component of the quasistanding-wave guided field formed in the cavity. The rates 
are normalized to the free-space decay rate 70 = u^d 2 / (3nheoc 3 ) . Figure [2]( a) shows that 
the cavity-modified rate of spontaneous emission into guided modes 7 cav g y d and the total 
spontaneous emission rate V increase with increasing reflectivity |-R| 2 . Meanwhile, the rate 
of spontaneous emission into radiation modes 7 ra d does not depend on |i?| 2 . In the absence of 
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the cavity (\R\ = 0), the rates of spontaneous emission into guided and radiation modes are 
7g y d = 7cav gy d(|-R| = 0) = O.4870 and 7rad — l-257o, respectively, and the total spontaneous 
emission rate is 7 = r(|i?| = 0) = 1.737 . Figure Mjoi) shows that the rate of spontaneous 
emission into guided modes 7 cav gyd and the total spontaneous emission rate F have a peak 
when the fiber radius a is around 191 nm. 




a 
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FIG. 3. Dependences of the channeling efficiency r] = 7 cav gyd/r on (a) the FBG mirror reflectivity 
|i?| 2 and (b) the fiber radius a for the parameters of Fig. [2J 

The efficiency of channeling of emission into guided modes is characterized by the param- 
eter 77 = 7 C avgyd/r. In Fig. Owe plot 77 as functions of the FBG mirror reflectivity \R\ 2 and 
the fiber radius a for the parameters of Fig. [21 Figure [3]^a) shows that the channeling effi- 
ciency 77 increases with increasing reflectivity |i?| 2 and can achieve substantial values when 
|i?| 2 is close to unity. Indeed, for the reflectivity |i?| 2 = 0.8 or 0.9, we obtain 77 = 0.87 (i.e. 
87%) or 0.94 (i.e. 94%), respectively. Figure [3]^b) shows that the channeling efficiency 77 
achieves a peak when the fiber radius a is around 191 nm. It is interesting to note that, due 
to the FBG cavity, the channeling efficiency 77 can achieve substantial values in a relatively 
wide range of a. Indeed, for |i?| 2 = 0.9, we find 77 > 80% when a is in the range from 130 to 
300 nm. 

The total spontaneous emission rate T and its components 7 ca vg y d and 7 ra d depend on 
the axial coordinate z and the radial coordinate r of the atom. In Fig. H]we plot T, 7 ca vg y d, 
and 7 ra( j as functions of z and r in the case where the dipole of the atom is oriented along 
one of the spherical basis vectors u±j.. Figure IU(a) shows that the cavity- modified rate of 
spontaneous emission into guided modes 7 ca vg y d and the total spontaneous emission rate T 
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(a) (b) 




FIG. 4. Dependences of the total spontaneous emission rate T (solid lines) and its two components 
Tcavgyd (dashed lines) and 7 ra d (dotted lines) on (a) the axial coordinate z and (b) the radial 
coordinate r of the atom in the case where the dipole of the atom is oriented along one of the 
spherical basis vectors u±i. In (a), the atom is located at r = a (on the fiber surface). In (b), 
the atom is located at z = (at the center of the cavity). The fiber radius is a = 200 nm. The 
reflectivity of the FBG mirrors is |i?| 2 = 0.9. Other parameters are as in Fig. [2j 

oscillate with varying z. The spatial period of the oscillations is tt/Po. The maxima and 
minima of the rate 7 ca vg y d or T correspond to the enhancement and inhibition, respectively, 
caused by the FBG cavity, and are achieved at the antinodes and nodes, respectively, of the 
parallel-to-dipole component of the quasistanding-wave guided field formed in the cavity. 
Meanwhile, the rate of spontaneous emission into radiation modes 7 ra( j does not depend on 
z and is finite. This explains the observation in Fig. H](a) that, at the nodes of the cavity 
field, the total atomic decay rate T remains finite although the component 7 ca vg y d becomes 
very small. Figure H](b) shows that the effect of the fiber on T, 7 cav g y d, and 7 ra d is largest for 
the atom on the fiber surface. It is clear that, when the atom is located at an antinode of 
the parallel-to-dipole component of the cavity guided field and is near to the fiber surface, 
7cav gy d and consequently T are substantially enhanced by the cavity. When the atom is far 
away from the fiber (r ^> a), the rate 7 cavg yd reduces to zero while the rates 7 ra d and V 
approach the free-space value 70. 

We plot in Fig. [5] the efficiency of channeling of emission into guided modes r\ = 7 cavg yd/r 
against the axial coordinate z and the radial coordinate r of the atom for the case of Fig. 
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(a) , (b) 




Po z r - a [run] 

FIG. 5. Dependences of the channeling efficiency r\ = 7 ca vgyd/r on (a) the axial coordinate z and 
(b) the radial coordinate r of the atom for the parameters of Fig. [H 

HI Figure [5]^a) shows that r\ oscillates with varying z, with the period 7t/Pq. It is interesting 
to note that 77 is substantial in broad regions around the antinodes of the parallel-to-dipole 
component of the cavity guided field and has narrow dips at the nodes. The appearance of 
such features is due to the fact that the total atomic decay rate T has two components, one 
is enhanced or inhibited around the antinodes or nodes, respectively, and the other is not 
modified by the cavity and is substantial. Figure EJ^b) shows that the channeling efficiency 
T) reduces with increasing atom-to-surface distance r — a and is substantial in a broad region 
of r — a. Indeed, more than 50% of emission can be directed into guided modes when the 
atom-to-surface distance is less than 350 nm. In addition, the channeling efficiency rj can 
be significant even when r — a is large. Indeed, up to about 15% of emission can be directed 
into guided modes when the atom-to-surface distance is 600 nm. 

The numerical results presented in Figs. EHSwere obtained for the case where the dipole 
vector d of the atom is perpendicular to the fiber axis z. Meanwhile, the spontaneous decay 
rates and their modifications caused by the FBG cavity depend on the orientation of the 
atomic dipole. We plot in Figs. [6] and [7] the spatial dependences of the decay rates T, 7 ca v g yd, 
and 7 ra d and the channeling efficiency 77 for the case where the dipole of the atom is oriented 
along the fiber axis z. Comparison between Figs. H] and [6] and between Figs. [5] and [7J shows 
that the decay rates and the channeling efficiency are smaller for an atom with a dipole 
parallel to the fiber axis than for an atom with a dipole perpendicular to the fiber axis. In 
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FIG. 6. Dependences of the total spontaneous emission rate T (solid lines) and its two components 
Tcavgyd (dashed lines) and 7 ra d (dotted lines) on (a) the axial coordinate z and (b) the radial 
coordinate r of the atom in the case where the dipole of the atom is oriented along the fiber axis 
z. In (a), the atom is located at r = a (on the fiber surface). In (b), the atom is located at 
(3qz = ±7r/2 (one- fourth of the guided- light wavelength from the cavity center). The fiber radius is 
a = 200 nm. The reflectivity of the FBG mirrors is \R\ 2 = 0.9. Other parameters are as in Fig. [2l 




200 400 600 

r - a [nm] 



FIG. 7. Dependences of the channeling efficiency r] = 7 ca vgyd/r on (a) the axial coordinate z and 
(b) the radial coordinate r of the atom for the parameters of Fig. [6j 

addition, we observe that the positions of maxima (minima) in the case of Figs. ID^a) and 
[5]^a) correspond to the positions of minima (maxima) in the case of Figs. [6]^a) and 0(a). 
Such opposite behaviors are due to the differences between the phase shifts per reflection of 
the longitudinal (q = 0) and the transverse (q = ±1) components of the guided field. 
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V. DELAY-DIFFERENTIAL EQUATION FOR MULTIPLE REFLECTIONS 



We now examine Eq. (fT2j) in the case where the coupling between the atom and the cavity- 
modified guided modes may be strong and, consequently, the Born-Markov approximation 
for the contribution of the cavity guided modes to the atomic decay may not be valid. 



We follow the approach of Refs. [4l|, |42j, |43| and derive a delay-differential equation that 



describes explicitly multiple reflections in our model. 

We start from Eq. (fT8|) for the cavity impact factor G(u). We expand the denominator 
of the fraction in this equation into a Fourier series as 

1 

1 - \R\ 2 + A\R\ 2 (1 - [Rl 2 )- 1 sin 2 $ 
2 ^ \ R 



\R\ 2n 

^T-^— cos(2n$). (29) 



l + \R\^ l + S n , 

With the help of the above formula, we expand the cavity impact factor G(uj) into a series 

as 

\_R\_ 

t i — ' I ^ It,. ' ! . . 

n=0 ' n=0 



\2n 

\2n+l 



G ^ = 2 E tVa — co < 2n ^ + E w 

„ =0 1 + n=0 

x { cos[(2n + 1)$ + 2/3z] + cos[(2n + 1)$ - 2/3z]}. 

(30) 

We insert Eq. ( 130]) into Eq. f|T6|) and calculate the integrals with the help of the formulas 

rOO 

cos(2n$)e~ i(< ^ o)T Aj 

-oo 

= (l + 5 nfi )7ie 2m ^5(T-2nr L ), 
cos[(2n + 1)$ ± 2[3z}e- i{uJ - UJo)T dLU 

-oo 

= 7re (2n+1)i *°e ±2i/3o ^(r - 2nr L - r±), (31) 

where r + = (L + 2z)/v g and r_ = (L — 2z)/v g are the position-dependent group delays 
due to the left and right mirrors, respectively, and t l = L/v g = (r + + r_)/2 is the group 
delay per cavity crossing. In deriving expressions fl3Tl) we have neglected the group velocity 
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dispersion. When insert the result of the calculations into Eq. fjl2j) . we obtain 



oo 



c a (t) 



7gyd 

2 



c a we(t) + 2^| J R| 2 V ni *° 



n=l 



x C (t - 2nr L )6(t - 2nr L ) + ^ | jR |2n+i e (2n+i)i* 



n=0 



x e 2i/3 °*C a (t - 2nr L - r+)6(t - 2nr L - r+) 



+ e - 2ip ° z C a {t - 2nr L - r_)6(t - 2nr L - r_) 




(32) 



Here B(i) is the Heaviside step function, whose value is zero for negative argument and one 
for positive argument. 

Equation ( 132]) is a delay-differential equation for the decay of an atom near a fiber with a 
pair of FBG mirrors. The first term, C a (t)Q(t), does not depend on the reflection coefficient 
R. This term describes spontaneous emission into guided modes in the absence of the cavity. 
The other terms are associated with the coefficients of the type R n , where n = 1,2,.... 
Such terms describe the backaction of the emitted photon on the atom after the photon is 
reflected from the mirrors n times. The quantities of the type 2nr^ and 2nr^ + r± are the 
group delays. The factors of the type e 2m ®° and e ( 2n + 1 ) l *o e ±2«/3 ^ describe the phase shifts 
of the parallel-to-dipole component of the guided field due to the propagation along the 
nanofiber and the reflection from the FBG gratings. Thus, the delay-differential equation 
( 132]) describes spontaneous emission of the atom in terms of multiple reflections. Due to 
the effect of retardation on the atomic state, the atomic decay may become nonexponential. 
We note that the absorption of the guided field by the fiber material can be incorporated 
into the theory by adding an imaginary part to the longitudinal wave number /3q, which 
appears in expression (123]) for the phase shift per cavity crossing $ and also in the local 
phase factors e *° z . 

It is clear from Eq. (132]) that, when t < r min = min{T£, r + , r_}, we have C a = —( f y/2)C a . 
Here, 7 = 7gyd + 7 ra d is the total rate of spontaneous emission into both types of modes in 
the absence of the cavity. The above result means that the atom does not feel the presence 
of the cavity until the time t = r m i n . 



In the framework of the Born-Markov approximation, we can replace the variables C a (t — 
2nr L ) and C a {t-2nT L -r±) in Eq. ® by C a (t). Then, Eq. ([32]) reduces to C a = -{T/2)C a . 
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Here, T = 7 gy dG ? o + 7 ra d is the total rate of spontaneous emission into both types of modes 
in the presence of the cavity. This result is in agreement with the results of Sec. IIVI on the 
exponential decay of the atom in the overdamped- cavity regime [see Eqs. (1241) and (I25|) ]. 

The delay-differential equation (13"!2"1) is similar to but different from the corresponding 
equation for the case of planar Fabry-Perot microcavities [43J . The key difference is that the 
coefficients in the delay-differential equation for planar Fabry-Perot microcavities contain 
£ _:L -, £~ 2 -, and £ -3 -terms which correspond to the dipole radiation field, the induced field, 
and the electrostatic field, respectively, due to the mirror images [43|. Here, £ = 2nu r L , 
2nujQTL + ojqt +) or 2uu tl + uj t^ is the retardation time. The absence of the £~ 2 -, 
and £ -3 -terms in the expressions for the coefficients in Eq. (1321) is because the FBG cavity 
reflects only the guided modes and is therefore similar to one- dimensional cavities. Due to 
this reason, Eq . (132]) is almost the same as the corresponding equation for one-dimensional 
cavities [4l|, |42j. A difference between the two cases is that Eq. (1321) contains an additional 
term, namely the term — (7 ra d/2)C a (t), which describes spontaneous emission from the atom 
into radiation modes. Another difference is that the cross-section area of the cavity modes 
is rigorously included in the expression for the rate 7 gy d of spontaneous emission into guided 



modes but is phenomenologically included in the treatment of Ref. 
treatment of Ref. 421 ] . 



411 ] or is omitted in the 



The delay-differential equation (132]) fo r the probability amplitude C a of the atomic upper 
state | a) can be solved numerically^, 1^2 ] by using a subroutine developed in Ref. 48]. 
We solve this equation and plot in Fig. El the time evolution of the atomic upper-state 
population P a = \C a \ 2 for the case where the FBG cavity length is L = 20 cm. The atom is 
located on the fiber surface (r = a) and at the cavity center (z = 0). The cavity length is 
tuned to resonance with the atomic transition frequency so that the phase shift per cavity 
crossing $ is an even or odd multiple of it, that is, the center of the cavity corresponds to 
an antinode or a node, respectively, of the parallel-to-dipole component of the cavity guided 
field. Other parameters are as in Fig. [21 For comparison, the exponential decay of the 
atomic upper-state population in the absence of the cavity is shown by the dashed curves. 



The solid line in Fig. 0(a) shows the occurrence of vacuum Rabi oscillations 4jJ, |42j, |43 |. 
Such oscillations are due to strong coupling between the atom and the guided field in the 
FBG cavity. It is interesting to note that strong coupling and vacuum Rabi oscillations 
can occur even when the cavity length is large (L = 20 cm) and the finesse of the cavity 
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FIG. 8. Time evolution of the upper-state population P a = \C a \ 2 of the atom in a long FBG 
cavity. The length of the cavity is L = 20 cm and is tuned to resonance with the atomic transition 
frequency so that the phase shift per cavity crossing <3?o is an (a) even or (b) odd multiple of ir. 
The reflectivity of the FBG mirrors is \R\ 2 = 0.9. The atom is located on the fiber surface and at 
the cavity center. Other parameters are as in Fig. [2J For comparison, the exponential decay of the 
atomic upper-state population in the absence of the cavity is shown by the dashed curves. 

is moderate [F = 7r|i2|/(l — |-R| 2 ) — 30]. There are two reasons for this. The first reason 
is that the field in the guided modes of the nanoflber is confined in a small area of the 
transverse plane, that is, the guided-mode cross-section area is small. Due to this reason, 
the effective cavity-mode volume can be small and, consequently, the cavity-atom coupling 
constant can be large even when the FBG cavity length is large. The other reason is that the 
FBG cavity is similar to one- dimensional cavities. In such a cavity, the cavity damping rate 
reduces with increasing cavity length faster than the strength of the coupling between the 
atom and the cavity. Unlike one-dimensional cavities, planar Fabry- Perot optical cavities 
have off-axis modes, which reduce the cavity QED effects 
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43j. In addition, the radius 



of the cavity mode in a planar Fabry-Perot optical cavity increases with increasing cavity 



length L and with increasing mirror reflectivity |i?| 2 40j. In the case of curved Fabry- 
Perot optical cavities, the typical values of the mode waist are much larger the wavelength 
of light. Consequently, the realization of strong coupling in a planar or a curved Fabry- 
Perot cavity requires a smaller cavity length and a higher finesse than in a nanofiber-based 
cavity. The typical Fabry-Perot optical cavities used in experimental realizations of strong 
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couplin g h ave lengths in the range from 10 to 100 /im and finesse factors on the order of 10 E 



3, E 



18 



19 



20 



2l|. 



The solid line in Fig. [3(b) shows the decay of the atom is almost exponential. Comparison 
between the solid and the dashed lines shows that the atomic decay is slightly slowed down 
by presence of the cavity. In the case of this figure, the atom is positioned at a node of the 
cavity guided field and, therefore, spontaneous emission of the atom into guided modes is 
substantially inhibited. The total spontaneous emission of the atom is mainly determined 
by spontaneous emission into radiation modes. This decay channel is slightly weaker than 
the cavity-free atomic decay. 
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FIG. 9. Time evolution of the upper-state population P a = \C a \ 2 of the atom in a short FBG 
cavity. The length of the cavity is L = 2 mm and is tuned to resonance with the atomic transition 
frequency so that the phase shift per cavity crossing $0 is an ( a ) even or (b) odd multiple of ir. 
The reflectivity of the FBG mirrors is \R\ 2 = 0.9. The atom is located on the fiber surface and 
at the cavity center. Other parameters are as in Figs. [2] and For comparison, the exponential 
decay of the atomic upper-state population in the absence of the cavity is shown by the dashed 
curves. 

We plot in Fig. [9] the time evolution of the atomic upper-state population P a = \C a \ 2 for 
the case where the FBG cavity length is L = 2 mm. Other parameters are as in Figs. [2] 
and[SJ The figure shows that the atomic population decay is almost exponential. According 
to Fig. [21(a), the exponential decay of the atom at an antinode of the cavity guided field 
(see the solid curve) is substantially faster than the cavity-free atomic decay (see the dashed 
curve). According to Fig. [9](b), the exponential decay of the atom at a node of the cavity 
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guided field (see the solid curve) is slightly slower than the cavity-free atomic decay (see the 
dashed curve). 

Comparison between Figs. 0(a) and 0(a) shows that vacuum Rabi oscillations can be 
observed in the case of Fig. Ela), where the cavity length is rather large (L = 20 cm), but 
not in the case of Fig. M^&), where the cavity length is much shorter (L = 2 mm). Thus, 
vacuum Rabi oscillations cannot occur when the cavity length is too short. This result is 
different from the common belief that the smaller cavity can produce the stronger vacuum 
Rabi oscillations [kJ. Such a belief was based on the results for high- finesse microcavities. 
Meanwhile, our model involves the use of a moderate-finesse nanofiber-based cavity. 

Comparison between Figs. E(b) and [21(b) shows that the time dependences of the atomic 
upper-state population P a in the two cases are essentially the same. Moreover, they are 
almost identical to the exponential decay of the atom into radiation modes. The reason 
is the following: in the two cases, the atom is positioned at a node of the cavity guided 
field and, hence, spontaneous emission into guided modes is inhibited. Since this effect is 
substantial enough, the total atomic decay process is mainly determined by the process of 
spontaneous emission into radiation modes. 

We plot in Fig. [10] the time evolution of the atomic upper-state population P a = \C a \ 2 
for different values of the cavity length, in the range from 100 m to 1 mm. The length 
of the cavity is tuned to exact resonance with the atomic transition frequency so that the 
phase shift per cavity crossing <3> is an even multiple of ir, that is, the center of the cavity 
corresponds to an antinode of the parallel-to-dipole component of the cavity guided field. 
Since the cavity length is rather large in the cases of Figs. fTOT a)-(c). we take into account the 
absorption of the guided light by the fiber material in the calculations. For this purpose, we 
add an imaginary part of a/2, with a = 10~ 5 cm -1 being the typical absorption coefficient 
for silica, to the longitudinal wave number /?o, which appears in Eq. (|32|) through the phase 
shifts $ an d ±2/3 z. However, we neglect the nonradiative atomic decay caused by the 



material absorption 49 



50] . The ratio of the rate 7nonrad of such a nonradiative process to 



the natural decay rate 70 is given in the limit of small atom-to-surface distances r — a by 
the factor e 7 /[2|e + l| 2 ^ (r — a) ], where ej is the imaginary part of the complex permitivity 



49 



501 ] . In the case of silica, ej is on the order of 10 10 . Therefore, the nonradiative 



decay rate 7nonrad of an atom with the transition wavelength Ao = 852 nm of the cesium D2 
line is significant only when the distance r — a from the atom to the fiber surface is on the 
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FIG. 10. Time evolution of the upper-state population P a = \C a \ 2 of the atom in a FBG cavity 
with length L = 100 m (a), 10 m (b), 1 m (c), 10 cm (d), 1 cm (e), and 1 mm (f). The length of 
the cavity is tuned to exact resonance with the atomic transition frequency so that the phase shift 
per cavity crossing <3?o is an even multiple of ir. The atom is located on the fiber surface and at 
the cavity center. The reflectivity of the FBG mirrors is \R\ 2 = 0.9. Other parameters are as in 

Fig. m 

order of or less than 0.2 A. Such a threshold distance is very small as compared to the light 
wavelength and the fiber radius, and is even smaller than the Bohr radius. Therefore, it is 
neglected in our treatment. The aim of the choice of the value r — a = for the calculations 
of Fig. [10] as well as Figs. El and [9] is to show the most dramatic effects in the limiting case 
where the effects of the material absorption, the surface-induced potential, and the surface 
roughness on the atomic decay can be neglected. Our additional calculations, not shown 
here, confirm that the numerical results presented in Figs. IBHTOl remain basically unchanged 
when the value zero for r — a is replaced by a few nanometers. 

Figure [TOT a) shows that, when the cavity length L is large enough, the decay of the atomic 
upper-state population P a is almost exponential. Such a decay is close to the exponential 
decay of the atom in the absence of the cavity. The cavity-free atomic decay rate is 7 = 
7gyd + 7rad and is approximately equal to 1.737o in the case considered. We note that, in the 
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case of large L, we may also observe vacuum Rabi oscillations. However, such oscillations 
are weak, and the corresponding period is large and approaches the cavity crossing time 7£. 

Figures [TuYc) and [TDlfd) show that significant vacuum Rabi oscillations occur when the 
cavity length L is in a range on the order of 10 cm to 1 m. Such lengths are rather large. 



We mention again that, in the case of planar 
18 
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431 ] and curved 
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2l| Fabry-Perot cavities, due to the substantial magnitudes of the cavity-mode 



cross-section areas and the effects of the off-axis modes, strong coupling cannot be realized in 
long cavities. Comparison between Fig. ITOT c) and fTOT d) shows that a decrease in the cavity 
length leads to a decrease in the vacuum Rabi period. This feature is in agreement with the 
fact that the vacuum Rabi frequency Q is proportional to the factor 1/y/rE = ^v g /L [see 
Eq. PHI) ], which characterizes the cavity mode density or the inverse of the cavity mode 
volume. 

Figure fTDT f) shows that, when the cavity length L is small enough, the decay of the atomic 
upper-state population P a returns to the exponential-decay regime, with a cavity- modified 
decay rate T = 7 gy dG'o + 7 ra d- We find T = 19.337 = II.2O7 in the case of the figure. 
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FIG. 11. Time evolution of the upper-state population P a = \C a \ 2 of the atom at different distances 
r — a = (solid line), 50 nm (dashed line), and 100 nm (dotted line) from the fiber surface in a 
FBG cavity. The cavity length is L = 10 cm. The atom is located at the center of the cavity. 
Other parameters are as in Figs. [2land[T0l 

The dependence of the time evolution of the atomic upper-state population P a = \C a \ 2 
on the distance r — a from the atom to the fiber surface is illustrated in Fig. [Hj The figure 
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shows clearly that vacuum Rabi oscillations can be observed even when the distance r — a is 
as large as 100 nm. The strong coupling between such a distant atom and the guided field 
is due to the effect of the FBG cavity. 

We conclude this section by presenting an analytical solution to the delay-differential 
equation (I3"2"j) in a particular case where the atom is at the center of the cavity, i.e., z = 0. 
In this case, Eq. ( 132|) reduces to 

C a (t) = - 7gyd W^Cait - nr L )Q(t - nr L ) 

n=l 

- \C a {t). (33) 

The above equation has been solved analytically in Ref. 43j. When we follow the result of 
Ref. y, we find 

oo 

C a (t) = e" 7i/2 \R\ n e ni *°e niTL,2 Q{t - nr L ) 

n=0 

Ei n P (t-nr L ) p 

(-^ fejfcj I ' (34) 



""1 )■ ■ ■ s^n 

where the inner sum is over all non-negative integers k%, fe, • • • , k n that satisfy the condition 
k\ + 2/^2 + ■ • - + nk n = n, and p = k\ + ki + ■ ■ • + k n is their sum. Expression (134"1) describes the 
time dependence of the probability amplitude C a of the atomic upper state \a). The results 
of calculations of expression fl34l) are in complete agreement with the numerical solutions of 
Eq. (J33D. 



VI. SINGLE-MODE CAVITY 



In order to get insight into our model, we approximate the delay-differential equation fl32|) 



under 
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;he single-mode cavity condition. For this purpose, we follow the procedures of Refs. 
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43j. We consider a cavity mode, whose frequency ui c is determined by the resonance 



condition $(co> c ) = mir. Here m is an integer number. We introduce the parameter 

$(a; c ) - $(u ) rmrj-jo 

which characterizes the detuning of the cavity mode frequency uj c from the atomic transition 
frequency u . It is clear that the separation between the cavity-mode frequencies u> c is 
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Au c = tt/tl- We rewrite the delay-differential equation ([32]) as 

s OO 

C a (t) = -l^hJ2e- MiA+K/2)TL C a (t - 2nr L ) 

^ n=l 



OO 



e -(2n+l)(iA+ K /2)r L 



x e(t-2nr L ) + (-l) m J2 

n=0 

W°* Ca ( t _ 2nr L - r + )0(t - 2nr L - t+) 
+ e- 2if3oz C a {t - 2nr L - r_)6(t - 2nr L - r_) J 



where 



- (36) 

k = — - 1 In \R\\ (37) 

is the cavity damping rate. Since the frequency separation between the cavity modes is 
Auj c = tt/tl, the cavity finesse is approximately given by F = Auj c /k = vr/(2| In \R\\). 

Although Eq. ([56]) is valid for an arbitrary integer number m, we choose such an integer 
number m for which the mode frequency u c is closest to the atomic transition frequency ujq, 
that is, the detuning A is smallest. We consider the case where the condition |A|tl <C 1 
is satisfied. This condition means that |A| <C Au c , that is, the cavity-atom detuning A is 
much smaller than the cavity-mode frequency separation Auj c . In this case, the effect of the 
cavity mode with the frequency u c on spontaneous emission of the atom is dominant over 
that of other cavity modes. Furthermore, we assume that \R\ = 1, so we have ktl <C 1. In 
addition, we assume that t ^> Tl and jtl <C 1. Under the above conditions, we can replace 
the sums in Eq. (136|) by integrals and hence obtain 

C a (t) = cos 2 \(3 c z + mix j2) 
t 



x 



j e - {iA+K /w-t') Ca{t > )dt > _ 2c B (t), (38) 



o 



where (3 C = (3{uj c ) = j3 + A/v g . When we differentiate the above equation with respect to t, 
we find the second-order differential equation 

~Q 2 /.. k\j~ 



K + 7 

O a -h I llA -t- 

where 



C a = 0, (39) 



Q = 2^p^-\ cos(f3 c z + mvr/2) | (40) 
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is the cavity-atom coupling constant. Note that Eqs. (139 ) and (HOI) are in agreement with 

42] . It is clear from Eq. (j3Dj) that 
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the results for high-finesse one-dimensional cavities 
the cavity-atom coupling constant ft is inversely proportional to the factor ^frl = ^L/v g , 
which effectively characterizes the cavity mode volume or the inverse of the cavity mode 
density. In addition, ft is proportional to the rate of spontaneous emission into guided 
modes v /7 gy d- Since the field in guided modes is tightly confined in the transverse plane, 
that is, the guided-mode cross-section area is small, the rate 7 gy d can be substantial when 
the atom is close to the fiber surface 4|. Therefore, ft can achieve substantial values even 
when the cavity length L is large. 

We analyze the case of exact cavity-atom resonance, where A = 0, that is, u c = u and, 
consequently, (3 C = (3q. In this case, Eq. (1591 reduces to 



C a + 



K + 7 



C a + 



ft 2 + «7 



C n = 0. 



(41) 



The initial conditions for the spontaneous emission process are C a {0) = 1 and C a (0) = —7/2. 
For these initial conditions, the solution to Eq. (j4~Tl) is found to be 



C a (t) 



where 



-(k+7)*/4 



A 



cosh(At/2) + ^— - sinh(At/2) 



(42) 



Below, we study several different regimes of the general solution (1421) . 

First, we consider the strong-coupling (underdamped-cavity) regime, where Q is suffi- 
ciently large that the condition 

2ft >k, 7 (44) 
is satisfied. In this regime, Eq. (1421) yields 



C a (t) = e- {K+ ^ t/4 cos(Qt/2). 
Hence, the population P a (t) = \C a (t)\ 2 of the atomic upper state \a) is found to be 



(45) 



P a (t) = e 



-(«+7)i/2 „„ c 2 



cos 2 (ftt/2). 



(46) 



The above solution descri 
coupling regime [4l|, |42j, |43| 



Des the occurrence of vacuum Rabi oscillations in the strong- 
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The strong-coupling condition (jUJ) can be rewritten as 

L 2 > L > L u (47) 

where 

L2 = lK2kd cos 2 (/3oZ + m7r/2)) 

7 

= ^ g In 2 \R\ (A8) 

1 4 7gyd cos 2 (/5 2 + m7r/2)' 1 ; 

Condition f j4"Tj) says that the strong-coupling regime can be realized only if the cavity length 

L is sufficiently small as compared to L 2 and is sufficiently large as compared to L\. It is 

clear that condition (1471) can be realized only if 

L 2 > Li. (49) 

When the atom is positioned at a node of the parallel-to-dipole component of the cavity 
guided field, we have cos(/?o-2 + mir/2) = 0, which leads to L 2 = 0, L\ = oo, and Q = 0. In 
this case, condition (|49j) and the strong-coupling condition (jUJ) cannot be satisfied. 

When the atom is positioned at an antinode, we have cos(/3oz + wtt/2) = ±1, which leads 

to 

q = 2 J^ (50) 



and 

16v fl 7 gyd 



2 



7 2 



^ = -^1^^!= . (51) 



In this case, condition (j4"9~|) can be rewritten as 



7gyd 



> J — !— ^ = . (52) 

7 8 16F 1 ; 

Condition fl52l) can be satisfied if the cavity-free channeling efficiency factor 7gy d /7 = 
7gyd/ (7gyd + 7rad) and the cavity finesse F are sufficiently substantial. Under condition 
( I52p . we can choose an appropriate cavity length L that satisfies condition (|4"9"1) for the 
strong-coupling regime at an antinode. 

Since the field in guided modes of the nanofiber is confined in a small area of the transverse 
plane, that is, the guided-mode cross-section area is small, the channeling efficiency factor 
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7 gy d/7 can achieve substantial values when the atom is close to the fiber surface J4]. In 
this case, condition (1521 can be satisfied for moderate values of the finesse F of the cavity. 
Furthermore, the cavity-atom coupling constant Q and the upper limit value L 2 can be 
large. Consequently, the strong-coupling condition (J33J) and its equivalent form (14T1) can be 
satisfied for large values of the cavity length L. It is interesting to note that, in our model, 
the cavity damping rate k, given by Eq. (|37|) . decreases faster with increasing L than the 
cavity-atom coupling constant Q, given by Eq. (150]) . Due to this fact, the upper limit L 2 
for condition (14"T|) is determined by the requirement 2Q 3> 7 but not by the requirement 



2Q S> k. This is a common feature of one-dimensional cavities 
that, in the case of planar 
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431 ] and curved 
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421 ] . We emphasize again 
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2l| Fabry- 



Perot cavities, due to the substantial magnitudes of the cavity mode cross-section areas and 
the effects of the off-axis modes, strong coupling cannot be realized in large cavities. The 
typical lengths of Fabry-Perot optical cavities used in experimental realizations of strong 
coupling are in the range from 10 to 100 /im If], 17, 18, 19, 20, 21 1. Such short cavities 
must have high finesse in order to achieve the strong-coupling regime. 

Unlike the upper limit value L 2 for the strong-coupling condition (H7J), the lower limit 
value Li is determined by the requirement 2Q ^> k. In the case of high-finesse cavities, 
where F ^> 1, we have L\ — > 0. However, when F is moderate, L\ can become large. Thus, 
strong coupling cannot be realized in a FBG cavity with a moderate finesse F if the cavity 
length L is too short. In such a cavity, the cavity damping rate n is much larger than the 
cavity-atom coupling constant Q. 

We discuss the possibilities of strong coupling and consequential vacuum Rabi oscillations 
in the cases of Figs. [8]^a) and [10J In the cases of these figures, the atom is positioned at 
an antinode of the cavity guided field, and the mirror reflectivity is \R\ 2 = 0.9 (the cavity 
finesse is F = 30). For the parameters of these figures, we find the critical values L 2 = 17 
m and L\ = 1 cm. It is clear that condition (T47j) is satisfied in the case of Fig. [8](a), 
where L = 20 cm, and in the cases of Figs. [TOTc) and fTOTd). where the cavity length is 
L = 1 m and 10 cm, respectively. This explains why vacuum Rabi oscillations are observed 
in the above-mentioned figures. Furthermore, we obtain (Q, k, 7)/7o — (7.97,3.51,1.73) 
in the case of Fig. [H^a), (Q, k,j)/j = (3.56,0.70,1.73) in the case of Fig. [TDT c). and 
(f2, k, 7) / 70 — (11.27,7.02,1.73) in the case of Fig. [TOTd). These parameters satisfy the 
strong-coupling condition (jSJ). For the free-space atomic decay rate 70 = 5.2 MHz of the 
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cesium D2 line, the cavity-atom coupling constant achieves the values £7 = 42, 19, and 
59 MHz in the cases of Figs. E](a), flOlfc). and [TOl d). respectively. Such values of Q are 
comparable to the values obtained in the experiments on realization of strong coupling in 



high-finesse Fabry-Perot optical microcavities 
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2l| 



Next, we consider the overdamped- cavity regime, where the condition 

/«>2fi,7 (53) 
is satisfied. In this regime, Eq. (1421) yields 

C a (t) = e- Tt ' 2 (54) 



and, hence, we obtain 



Here 



P a (t) = \C a (t)\ 2 ^e~ Lt . (55) 



T = Tgyd^o + 7rad (56) 

is the total decay rate of the atom, with 

Q 2 2 

G = 1 + = 1 + | r-r cos 2 faz + mvr/2) (57) 

K7gyd |l n |-R|| 

being the cavity impact factor for the rate of spontaneous emission into guided modes. The 
maximal enhancement factor is 

2 AF 

G ma x = 1 + tt—^t = 1 + — (58) 

| in | | | 7T 

and the maximal inhibition factor is 

G min = 1. (59) 

We note that Eqs. (EZJ) ([59]) agree with Eqs. (ISBVpBl of Sec. IV]in the limit \R\ 1. 
The overdamped-cavity condition ( |53l) can be rewritten as 



L< J L 1 ,L 3 , (60) 

where 

£3 = ^ Nfl il^. (61) 

7 1 1 ^7 
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Condition (I60j) indicates that the overdamped- cavity regime, where the spontaneous emission 
of the atom is an exponential-decay process with a cavity-modified rate T, can be realized 
only when the FBG cavity is sufficiently short. When the finesse F of the cavity is moderate, 
Li and L 3 can be large. For the parameters of Figs. Oanddni we find L\ = 1 cm and L 3 = 41 
cm. Then, the overdamped-cavity condition (I53"|) becomes L C 1 cm. It is clear that the 
case of Fig. [9j where L = 2 mm, and the case of Fig. [TUlff). where L = 1 mm, correspond 
to the overdamped-cavity regime. 

Finally, we discuss the case where L ^> L 2 , L 3 . In this case, we have 7 ^> Q, k. Then, Eq. 
(I42p yields C a (t) = e -7 ^ 2 and, hence, we find P a (t) = e -7 *. Thus, when the cavity is very 
long, the upper-state population P a of the atom reduces exponentially with the cavity-free 
atomic decay rate 7. Such a decay is observed in Fig. HUT a) although the parameters for 
this figure do not satisfy the conditions t ^> tl and <C 1, which were used in deriving 
Eq. flSHD from Eq. (JHBD - 

VII. SUMMARY 

We have studied spontaneous emission of an atom near a nanofiber with two fiber Bragg 
grating (FBG) mirrors. We have shown that the coupling between the atom and the guided 
modes of the nanofiber can be significantly enhanced by the FBG cavity even when the cavity 
finesse is moderate. We have found that, when the fiber radius is 200 nm and the cavity 
finesse is about 30, up to 94% of spontaneous emission from the atom can be channeled into 
the guided modes in the overdamped-cavity regime. 

We have derived a delay-differential equation which explicitly describes the effects of 
multiple reflections of the guided field on the atom. We have analyzed this equation in 
different regimes of the atomic decay. We have shown numerically and analytically that 
vacuum Rabi oscillations and strong coupling can occur in the FBG cavity even when the 
cavity finesse is moderate (about 30) and the cavity length is large (on the order of 10 cm 
to 1 m), unlike the case of planar and curved Fabry-Perot cavities. We have identified two 
reasons for this possibility. One reason is that the field in the guided modes of the nanofiber 
is confined in a small area of the transverse plane. Due to this reason, the effective cavity- 
mode volume can be small even when the FBG cavity length is large. Another reason is that 
the FBG cavity is similar to one- dimensional cavities, where there are no off-axis modes. 
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Appendix A: Mode functions of the fundamental guided modes of a nanofiber 



For the fundamental guided modes, the propagation constant (5 is determined by the fiber 



eigenvalue equation 



37] 



J (ha) 
haJi(ha) 



n\ + n 2 2 K[(qa) 1 



2n\ qaKi(qa) h 2 a 2 



n\ K[(qa) 



+ 



/i 2 



2n\ qaKi(qa) 



1 1 

+ 



1 1/2 



(Al) 



nfk 2 \q 2 a 2 h 2 a 2 

Here the parameters h = [n\k 2 — /3 2 ) 1 / 2 and q = ((3 2 — n\k 2 ) 1 ^ 2 characterize the fields inside 
and outside the fiber, respectively. The notation J n and K n stand for the Bessel functions 
of the first kind and the modified Bessel functions of the second kind, respectively. 

The mode functions of the electric parts of the fundamental guided modes 3TJ] are given, 
for r < a, by 

flKi(ga), 



» 



h J 1 (ha] 



■[(l-s)J (hr) 



s)J 2 (hr)\, 



e M = _ZC?^M[(1 - s)J (hr) + (1 + s)J 2 (hr)], 



fC 



h Ji(ha) 
2g K l {qa) 
(3 Ji(ha) 



Ji(hr), 



(A2) 



and, for r > a, by 



» 



er = iC[{\ - s)K (qr) + (1 + s)K 2 (qr)\, 
ef = -IC[{\ - s)K (qr) - (1 + s)K 2 (qr)} 



fC 2 A Kl {qr). 



(A3) 



Here the parameter s is defined as s = {l/q 2 a 2 +l/h 2 a 2 )/[J[{ha)/haJi{ha)+K[{qa)/qaKi{qa)\. 
and the coefficient C is determined from the normalization condition 

»2tt 



dip 



n 2 { |e (M) | 2 r dr = 1. 



(A4) 



Here n v Ar) = rii for r < a, and n r Ar) = n 2 for r > a. 
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Appendix B: Mode functions of the radiation modes of a nanofiber 



For the radiation modes, we have — kri2 < (3 < kn2- The characteristic parameters for 
the field in the inside and outside of the fiber are h = \Jk 2 n\ — (3 2 and q = \Jk 2 n\ — (3 2 , 
resp ectively. The mode functions of the electric parts of the radiation modes v = (uflml) 
371 ] are given, for r < a, by 



h? 



» - _ 



v h 2 

= AJ m (hr) 



(3hAJ L(hr) + im—^-BJ m (hr) 
r 

im—AJ m {hr) — hu>fj, BJ'(hr) 
r 



(Bl) 



and, for r > a, by 



^ = i E W-flSV) + 
q 3=1,2 



im- 



/m 



■^H^iq^-qu^DjH^'iqr) 



4 3=1,2 
3=1,2 

The coefficients Cj and £)j are related to the coefficients A and B as 



An 2 . 



-lY'-^iALj + ifXocBVj), 

l 2 



Dj = (-1) 



ieocAVj - BMj), 



(B2) 



(B3) 



where 



V 3 = ^(n 2 2 -nl)J m (ha)H£>(qa), 



Ma 



h 



J' m (ha)H£h(q a ) - -J m (ha)H$*'(qa), 



Q 



4 = ^JL(ha)H£>*(qa) - ^-J m (ha)H«>\qa). 



(B4) 



We specify two polarizations by choosing B = irjA and B = —irjA for I = + and / 
respectively. The orthogonality of the modes requires 



2~ 



dip 



n 



rf 



e M e ("> 



JO 

NJ w S(uj - u'). 



rdr 



J p=/3',m=m' 



(B5) 
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2 + 


\Lj\ 


2 


\Vj 


\ 2 + n 2 2 \ 


M 3 


2 ' 



This leads to 

V = ^ ^[ l 2 l7i\, - (B6) 
The normalization constant N v is given by 

= ^ + f|i>/) ■ (B7) 
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